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Abstract 

In the experiment of the quantum mirage, confinement of surface states in 
an ehiptical corral has been used to project the Kondo effect from one focus 
where a magnetic impurity was placed, to the other empty focus. The signa- 
ture of the Kondo effect is seen as a Fano antiresonance in scanning tunneling 
spectroscopy. This experiment combines the many-body physics of the Kondo 
effect with the subtle effects of confinement. In this work we review the essen- 
tial physics of the quantum mirage experiment, and present new calculations 
involving other geometries and more than one impurity in the corral, which 
should be relevant for other experiments that are being made, and to discern 
the relative importance of the hybridization of the impurity with surface (Vg) 
and bulk (V5) states. The intensity of the mirage imposes a lower bound to 
Vg/Vb which we estimate. Our emphasis is on the main physical ingredients 
of the phenomenon and the many-body aspects, like the dependence of the 
observed differential conductance with geometry, which cannot be calculated 
with alternative one-body theories. The system is described with an Ander- 
son impurity model solved using complementary approaches: perturbation 
theory in the Coulomb repulsion U, slave bosons in mean field and exact 
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diagonalization plus embedding. 
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I. INTRODUCTION 



The study of many-body phenomena in nanoscale systems is attracting a lot of attention 
in recent years. Progress in nanotechnology made it possible to construct nanodevices such 
as quantum dots (QD's) which act as ideal one-impurity systems in which the Kondo physics 
is clearly displayed [1-4] . The spectral density of localized electrons of a magnetic impurity 
in a metallic host, described by the impurity Anderson model [5], is known to display a 
resonance near the Fermi energy in the localized (or Kondo) regime [6]. The conductance 
through a QD is proportional to this density and calculations of the latter using the Wilson 
renormalization group leads to a good agreement with experiment [7,8]. In contrast to the 
density of localized electrons, the density of conduction electrons coupled to the former shows 
a dip or Fano antiresonance [9] (see section 3) . Using the fact that the tip of the scanning 
tunneling microscope (STM) captures essentially conduction electrons, Fano line shapes have 
been observed using scanning tunneling spectroscopy (STS) for different cases of magnetic 
impurities on metal surfaces [10-19]. They should also manifest in the conductance through 
quantum wires side-coupled to QD's [20-25]. The impurity Anderson model also describes 
the physics of the conductance through arrays of QD's weakly coupled to conducting leads 
[26]. 

A quantum corral is an area of about 40 nm^ delimited by a closed line of typically several 
tenths of atoms or molecules placed next to each other one at a time on an atomically 
flat metallic surface using a STM. The same microscope can be used to perform STS to 
study with meV resolution the electronic density inside these corrals [12,27-29]. Particularly 
interesting are the (111) surfaces of Cu, Au and Ag. These metals have nearly spherical 
Fermi surfaces with eight necks at the [111] and equivalent directions in which a gap opens. 
This allows the presence of Shockley states locahzed at the (111) surface uncoupled to 
the bulk states for small wave vector parallel to the surface and with nearly free electron 
dispersion [30]. STS experiments performed on these surfaces reveal fascinating standing- 
wave patterns and one can see the density coming from the wave functions obtained solving 
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the Schrodinger equation for a two-dimensional free electron gas inside a hard-wall corral 
[12,13,28,29,31]. Experiments in which different atoms or molecules were used to build the 
corral (Co, Fe, S, CO) suggest that the details of the boundaries are not important for the 
physics. A continuous description of the boundary is justified by the fact that the Fermi wave 
length I'KjkF ~ 3 nm is larger than the distance between adatoms. However, as discussed 
in section 6, the corrals are leaky and the hard- wall assumption should be abandoned for a 
quantitative description. 

The experiment of the quantum mirage is a beautiful combination of the physics of 
the quantum corral and the many-body Kondo effect. One Co atom acting as a magnetic 
impurity is placed at the focus of an elliptical corral built on the Cu(lll) surface, and 
a Fano dip is observed not only at the place of the impurity, but remarkably also at the 
empty focus with reduced magnitude [12]. Variants of this experiment involving other corral 
shapes and more that one impurity were presented in a conference [13]. In the original 
experiment, the space dependence of t^dl jdV ^ the differential conductance after subtracting 
the corresponding result without impurity, clearly resembles the density of the state number 
42 in increasing order of energy of free electrons in a hard- wall elliptical corral. This suggests 
that the main features of this space dependence can be explained by a one-body calculation. 
In fact, important features, like the possibility of obtaining mirages out of the foci can 
be understood by a simple tight-binding calculation [32] or from Green functions using 
hard- wall eigenstates [33]. Interesting effects like anti mirages were predicted for a non- 
magnetic impurity inside a hard- wall elliptical corral [34], and quantum mirages in s-wave 
superconductors were calculated [35]. Also phenomenological scattering theories in which 
the energy dependence of the Kondo resonance (directly related with the voltage dependence 
of dl /dV) as well as an inelastic part of the scattering are taken from experiment, are able 
to describe quantitatively the space dependence [31,36,37]. However, the calculation of the 
line shape of dl/dV, its dependence on the particular geometry of the corral, temperature 
or magnetic field, and the effects of interaction between impurities is out of the scope of 
these one-body approaches. 

4 



The first calculation of the voltage dependence of dl/dV has been done by one of us 
using perturbation theory in the Coulomb repulsion of the Anderson model U [38,39] A 
many-body calculation of the mirage effect is a challenge due to the particular nature of 
the conduction states brought by the confinement in the corral. In particular, available 
exact results for thermodynamic properties of the Kondo and Anderson impurities, obtained 
with the Bethe ansatz, assume a constant density of conduction states [40-42], while the 
Wilson renormalization group [43] (which allows accurate calculation of dynamical properties 
and was used in the context of nanoscopic systems and STS [7,8,44-46]), requires high 
symmetry around the impurity. If only a finite number of hard-wall eigenstates with well 
defined energies are mixed with the impurity (a problem that can be treated with exact 
diagonalization [47,48]), the hne shape of AdI/dV becomes qualitatively wrong (see section 
6). The reason is that the separation of the relevant energy levels is large in comparison 
with the Kondo temperature Tk ~ 53 K, while one knows that for a well developed Kondo 
resonance to exist, should be larger than the average separation of the relevant levels [49]. 
This points to the need of including a finite width of the corral eigenstates, which become 
resonances [38]. This need persists in presence of direct hybridization of the impurity with 
bulk states as shown in section 6 [50]. The width of the resonances cannot be too large 
because otherwise the space dependence of the state 42 of the elliptical corral, observed in 
dl/dV [12] would be blurred. 

A subject of current interest and debate in the literature is the relative importance of 
the hybridization of the impurity with bulk Vf, and surface Vs states. A first principles 
calculation seems not possible because of the large supercells needed. They should contain 
more than 10 layers perpendicular to the [111] direction in order for the Shockley surface 
state to develop, and more than 100 atoms per layer to reach the dilute limit of Co impurities 
on the surface [54]. On the basis of the rapid decay in AdI/dV as the STM tip is moved 
away from an impurity on a clean (111) surface, and a jellium theory of Plihal and Gadzuk 
[55], Knorr et al. concluded that bulk states dominate the formation of the Kondo singlet 
[17,18]. This is in agreement with tight- binding calculations [54]. However, recently Lin, 
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Castro Neto and Jones, using a nearly free electron approximation, including the effect of 
the gaps in the [111] and equivalent directions and calculating the wave functions under 
an adequate surface potential, concluded that the Kondo effect in the Cu(lll) surface is 
dominated by surface states [56]. They also obtained good agreement with experiment for 
the distance dependence of the amplitude of dl/dV and its voltage dependence on top of 
the impurity. Using a similar approach, but without attempting to solve the many-body 
problem. Merino and Gunnarsson concluded that surface states play an important role in the 
differential conductance for a system with a magnetic impurity on a clean (111) surface [57]. 
Therefore, the issue of the relative importance of V(, and Vg remains unclear. In contrast, in 
absence of the impurity, the relative contribution of the surface states to the conductance 
(STM tip-substrate hybridization) is known to be between 1/2 and 2/3 from experiments in 
which the bias voltage is swept below the bottom of the surface band (~ 0.45 eV below the 
Fermi energy) [17,58,59]. 

Since from the experiments we know that the presence of the corral strongly affects 
electronic structure of the surface states, it is clear that the variation of the line shape of 
dl/dV for different corrals or positions of the impurity inside the corral and its comparison 
with theory should help to elucidate the relative role of the hybridization of the impurity with 
surface and bulk states. A stronger sensitivity to the geometry implies a greater participation 
of the surface states in the formation of the Kondo resonance. Also the interaction between 
magnetic impurities inside a quantum corral should increase with the relative importance 
of surface states [47]. Unfortunately only the voltage dependence of dl/dV for a Co atom 
on a clean Cu(lll) surface and on an elliptical corral built on that surface is available 
for comparison [12]. Using perturbation theory in U both line shapes are qualitatively 
explained without bulk states [39] (see section 8). However, as we will show, this seems to 
be a particular case and usually the shape and width of the Fano dip are more sensitive to 
the geometry. 

In this work we discuss the main aspects of the physics of the quantum mirage. The 
emphasis is on the basic understanding of the phenomenon and its many-body aspects rather 
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than on quantitative fits. The latter would require more detailed knowledge of matrix 
elements and their wave vector dependence, crystal fields and other details. We extend 
previous many-body calculations for the space and voltage dependence of AdI/dV to new 
different situations. This can serve as a basis for comparison with experiment and help 
to elucidate the relative participation of surface and bulk states in the formation of the 
Kondo singlet for a Co atom on a Cu(lll) surface. We use three different many-body 
techniques: perturbation theory in U [60,61], exact diagonalization plus embedding [62-64] 
and a slave-boson mean-field approximation (SBMFA) [66,6,21,25,26]. The former two have 
already been applied to the quantum mirages [38,39,52,53,47] but have the disadvantage 
that they do not reproduce the correct exponential dependence of with the coupling 
constant for large [//A, where A is the resonance level width [65]. Therefore, the SBMFA 
is more appropriate to study the dependence of the width of the resonance on geometry. 

The paper is organized as follows. In section 2 we present the impurity Anderson model 
for either the corral or open surfaces, and discuss its assumptions and limitations. Section 3 
discusses the Kondo resonance and Fano antiresonance in the simplest version of the model 
for later comparison. The formalism and basic equations that determine the tunneling 
current are presented in section 4, using a many-body formalism, including tunneling of 
the tip of the STM with surface, bulk and impurity states. Section 5 is rather technical 
and explains the different many-body approaches. In section 6 we explain the effects of the 
confinement on the surface states, and how they are transmitted to the Kondo resonance 
and the line shape of the mirage effect. Section 7 is devoted to the space dependence of the 
differential conductance dl/dV inside an elliptical quantum corral, the effect of the impurity 
on it (Adl/dV), and the relation of these quantities with the wave functions of the surface 
states inside the corral. This brings insight into the effect of the width of the surface states 
and what controls the intensity at the mirage point. In section 8 we present results for 
the dependence of AdI/dV on bias voltage in different situations: clean surface, eUiptical 
corrals and a circular corral. In section 9 we estimate a lower bound for the participation of 
surface states in the Kondo resonance. In section 10 the interaction between two Anderson 
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impurities inside an elliptical corral is discussed. Section 11 contains a summary and a 
discussion. 



II. THE MODEL 



In this section, we explain and discuss the model used to describe the electronic structure 
of a system composed of one magnetic impurity interacting with surface and bulk states. 
The case of two impurities is left for section 10. The surface states can correspond to 
eigenstates of a clean perfect surface, or to a surface with a soft- wall corral. In both cases, 
the energy spectrum of the surface states is continuous. The wave functions fj{r) of the 
surface eigenstates are normalized in a large area [52,56]. Of course, all physical results are 
independent of this area. 

We take only one localized d orbital for the impurity. Technically it renders some many- 
body techniques easier (except the SBMFA). Previously Ujsaghy et al. [67] assumed a fully 
degenerate ground state while other recent calculations for impurities on (111) surfaces con- 
sidered the d322_j.2 orbital more important [56,57,68]. Tight binding calculations suggest 
that Cu(lll) surface states hybridize more strongly with the Co ?id^z^_r^ orbital, while bulk 
states prefer M^z and ?>dyz orbitals [54]. Recent accurate calculations using the Wilson 
renormalization group indicate that for the expected filling of near one d hole per Co im- 
purity, the Kondo resonance becomes strongly asymmetric in the orbitally degenerate case 
[46]. Then, one would expect in this case also a strongly asymmetric AdI/dV in contrast 
to the experimental observations for the (111) surface [12,17]. We neglect the s,p orbitals 
of the impurity. 

The Hamiltonian can be written as 



Sj, + e:)b\,bj, + -Bi ^ did, + Ud\d^d\di 



ja a 




(1) 



where s]^ (6]^) are creation operators for an electron in the j*'* surface (bulk) conduction 



8 



eigenstate in the absence of the impurity, but including the corral if present. The impurity 
is placed at the two-dimensional position Ri on the surface, and we assume that the hy- 
bridization of the impurity d orbital with the surface state j is proportional to its normalized 
wave function at that point (fij{Ri) [33,38]. Similarly for the bulk states, the hybridization 
is proportional to some average ipj{Ri) of the bulk wave function in the direction normal to 
the surface, that depends on 

Vi = VsX^^ (Ri), = V,i;j (Ri), (2) 

where Vg, H are energies representing local hybridizations in a tight- binding model [32,39] 
(see next section) and A = 2.38 A is the square root of the surface per Cu atom of a Cu(lll) 
surface. We assume also a constant density of bulk states. However, we must warn that 
recent calculations obtain a significant dependence of the matrix elements with wave vector 
[56,57,68]. This dependence affects the hne shape of dl/dV. Nevertheless, one expects that 
the trends in the modifications of the voltage dependence of dl/dV due to the modifications 
of the geometry remain the same, at least on a qualitative level. This is also suggested by 
the weak dependence of the results on the cutoff for the surface states Ec, which should be 
introduced in any theory of quantum corrals to avoid divergences in the Green's functions 
for the surface states. This can be though as an energy dependent hybridization Vg which is 
constant below Ef. and goes to zero abruptly at Ec- A hnearly decreasing Vg has also been 
used [38]. Increasing E^ leads to a weak increase in the width of the resonances and to a 
more asymmetric line shape, but the main conclusions regarding the mirage effect are not 
altered. 

The many-body part of the Hamiltonian which renders it non-trivial in all cases is the 
on-site Coulomb repulsion at the impurity site U. Another difficulty is the calculation of 
the surface wave functions (pj{r) for soft walls. They have been calculated exactly for a soft 
circular corral [52,53] and can be reasonably well approximated for an elliptical corral [52]. 
We return to this point in section 6. For open structures, the surface Green's function can 
be calculated using scattering theory [31]. However, this renders the many-body problem 
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too difficult. 



III. SIMPLE PICTURE OF THE KONDO RESONANCE AND FANO 

ANTIRESONANCE 

Before discussing the many-body techniques and the effect of the corral, we want to 
illustrate some basic features of the Kondo physics, using the simplest case of the Anderson 
model in which the impurity is hybridized with only one band (either surface or bulk) with 
constant density of states po ^-nd wave vector independent hybridization. The Hamiltonian is 
given by Eq. (1) eliminating the terms with Sj^, considering that the create Bloch waves, 
i^ka ~ (1/V^) S;exp(— i/c • Ri)bi^, where b]^ creates an electron at site / with position Ri) 
and taking — V/ y/N where N is the number of sites. Then, the impurity is hybridized 
with the band at one site that we call i. 

The impurity spectral density pd^iuj) of this model has been calculated accurately using 
Wilson renormalization group [69] and agree qualitatively with those of perturbation theory 
in U [61]. The results presented in Fig. 1 (a) were obtained using a self-consistent approach 
[70] based on an interpolation for the self energy of the Green's function between the ex- 
pression up to second order perturbation theory in the Coulomb repulsion U [60,61] and the 
exact result for [/ ^ oo. The resonant level width is A = npoV'^. This approximation works 
well for U <8A [23]. 
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FIG. 1. (a) spectral density of an impurity hybridized with a featureless conduction band and 
(b) spectral density at the conduction site hybridized with the impurity, as a function of energy at 
different temperatures. Parameters are — ep = — 3A and U = 8A. 

As seen in Fig. 1 (a), pdaiy^) shows characteristic charge fluctuation peaks (or shoulders 
for small U) at and + U and another peak near the Fermi energy ep characteristic 
of the Kondo regime. This peak is the so called Kondo resonance. Its half width at half 
maximum corresponds to the Kondo temperature T^-. At temperatures above the Kondo 
effect disappears and the spectral weight of the Kondo peak is transferred to the other two. 

The STS is much more sensitive to the conduction electrons than to the locahzed ones 
because the former are more extended in space and reach the tip of the STM with a larger 
amplitude. The site most affected by the impurity is the one which hybridizes with it (i). 
Using equations of motion (in the same way as done in the next section), it is easy to show 
that the Green function at this site is (the spin index is dropped for simplicity) 

Guioo) = GUuj) + [VGUu)YGa{uj), (3) 

where Gf^{uj) is the corresponding Green's function in the absence of the impurity and Gd{uj) 
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is the Green's function of the impurity. This equation is exact and does not depend on the 
approximations for G^di'-^)- If the unperturbed conduction band extends from —1^ to 1^ 
then 



W + CU 

\n\Z7—-\-t7rQ(W-U) 

W — LU 



(4) 



k 

If (as usual) W ^ and we are interested in energies |ci; — ei?| ~ T^, we can neglect the 
first term inside the brackets and approximate G°j(a;) = —iirpo. Replacing this in Eq. (3) 
and using Pia{oj) — — Im[G'ij(a;)]/7r one obtains the very simple result 

= 1 - npdaiuj)^- (5) 

Po 

Thus a peak in Pdai'^) implies a dip in the conduction density of states and more pronounced 
near the impurity (see Fig. 1 (b)). In more complex situations, in particular when confine- 
ment due to the corral is important, the real part of G%{cu) cannot be neglected and the dip 
in Piaifjj) is not directly related with the Kondo peak in pda{uj). In extreme cases, either the 
dip is replaced by a peak [56] or the structure near ep disappears, as we will show in section 
6. 

In any case, the above simple picture corresponds to a first rough approximation of the 
experimental observations of the voltage dependence of dl/dV for impurities on the (111) 
surfaces of Cu and noble metals, and we will use it for later comparison. 



IV. THE TUNNELING CONDUCTANCE 

In this section we write the basic equations which relate dl/dV with the Green's function 
at the impurity site. We include the hopping of the tip of the STM with the impurity, surface 
and bulk states in a many-body formalism. The tunneling geometry and energy diagram is 
shown for example in Fig. 1 of Ref. [31], but the impurity should be included if the tunneling 
current is measured near it [71], and also the bulk states according to experiment [17,58,59]. 

The total system S consists of a subsystem Sh described by the Hamiltonian H [Eq. (1)] 
and St contains the tip, which we assume can be described as a non-interacting system with 
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one-particle energies and the Fermi energy set at zero. Sh has all one particle energies, 
including the Fermi level, displaced by eV to lower energies by a bias voltage V, where e 
is the elementary charge. For simplicity we treat the case of positive V in which electrons 
are transferred from the tip to Sh- Extension to negative V is trivial using an electron-hole 
transformation. We assume a local hopping of the tip with the different states 

Hmia^^Aj^itLha + H.c), (6) 
fecr 

K ^\^Vj{Rt)sj^+pY,'^3iRt)hjc + q{\Rt - Ri\)da- (7) 

3 3 

Here creates an electron in the tip eigenstate k with spin cr, Rt describes the coordinates 
of the tip on the plane, and A, p, and q are parameters that describe the hopping of the 
tip with the different states of Sh- The function q{\Rt — Ri\) is small and decays strongly 
with the distance between the tip and the impurity \Rt — Ri\ due to the strongly localized 
nature of the impurity wave function. However, when i?t — i?j = 0, a small q introduces an 
important source of asymmetry in the line shape of dl/dV in addition to that corresponding 
to the structure of the Green's functions (see section 8). 

Treating Hmix in lowest order in perturbation theory and at T = 0, using Fermi's golden 
rule, the current due to the transfer of electrons from Sf to Sh becomes 

I = E \{^\J2^ltk.\9mE^ -E,- eV\ (8) 

V ka 

where \v)^ E^, are the eigenstates and energies of S and \g) is the ground state assumed non 
degenerate. Using the same notation with a subscript H for Sh-, we have \g) — rifeo- ^lo-lfl'-ff) 
with the product restricted to k such that < = 0. Replacing above and doing the 
calculations within St one has 

^ = E E \k^H\h\\9H)\^b{E^, - 6, - E,, - eV). (9) 

l^H fee 

Using the Lehman representation [72], the sum over vh is seen to represent the part of the 
spectral density pha{^k + eV^) of /ij. for electron addition. This corresponds to excitations 
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above ep with positive argument of pha [73]. By symmetry it is independent of cr. Trans- 
forming the sum over the tip states k as an integral assuming a constant density of states 
pt one gets 

/47re 

Il—A^ptj = j ^pha{e + eV)de = Pha{^)(^- (10) 

Prom here, it is clear that the differential conductance is proportional to the spectral density 
of the state ha{Rt)- 

dl/dV ^ PhaieV) = --ImGhaieV), (11) 

TT 

where Gha{^^) = {{ha] hl))u) is the Green's function of ha{Rt)- Therefore, in the rest of the 
paper we will be mainly concerned on the space and energy dependence of pha- This spectral 
density can be related with the Green's function for the d electrons Gdcr{oj) — {{da-; dl))^, 
and the unperturbed Green's functions for conduction electrons using equations of motion. 
Writing Cj^ to represent either sja- or bj^, the relevant equations can be written in the form 

{u;-e'^){{c,a;dl))^^Vj{{da;dl)U 

- e 4.))- = Vc{{da, 4)).- (12) 

Dropping the spin indices, using these equations and introducing the non-interacting Green's 
functions (in absence of the impurity) for conduction electrons 

ipj{Ri)ipAR2) 



G°(i?i,i?2,a;) = {{J2'PjiRi>j'r;¥AR2)s]j). = Yl 



GtiRM-E^-^^^^^ (13) 
■^r"^ — e,- 

the Green's function for the h operators becomes 

Gh{Rt:Ri:U;) = X^G^'.iRt^Ruu;) + p'^G^iRt, Rt,u;) + AGh{Ru Ri^u;), 
AGh{Rt, Ri, oo) = F{Rt, Ri, uj)F{Ri, Rt, uj)Gd{uj), 

F{R,,R2,u;)^VsX^G'i{Ri,R2,u;)+pVbGl{Ri,R2,u;) + q{\Ri-R2\). (14) 
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Here, the first two terms when replaced in Eq. (11) describe dl/dV in the absence of the 
impurity, while AG^ describes the effect of the impurity on the differential conductance 
Ml/dV. 

Note that the space dependence of dl/dV is determined only by the non-interacting 
conduction electron Green's functions. In particular at a distance of the impurity larger 
than ~0.5 nm, q{\Rt — Ri\) becomes irrelevant, the bulk part becomes less important in 
AdI/dV due to its more rapid decay with the distance between the tip and the impurity 
\Rt — Ri\ [74], and the space dependence is dominated by G%Rt, Ri,uj). The impurity 
Green's function Gd can only alter the relative weight of the real and imaginary part of the 
other factors in AG^. There is a natural length scale in the Kondo problem ^ = Hvf/Tk, 
where vp is the Fermi velocity. It has been interpreted as the size of the cloud of conduction 
electrons that screen the localized spin in the Kondo effect. The existence of this cloud 
is still controversial [75-77]. Theoretical work has shown that the persistent current as a 
function of flux in mesoscopic rings with quantum dots changes its shape smoothly as 
the length of the ring L goes through ^ and that jL is a universal function oi L/^ [78,79]. 
However, in our case, it is clear that ^ plays no role in the space dependence of dl/dV. 

V. THE MANY-BODY TECHNIQUES 

The core of the many-body problem is to solve the impurity Green's function Gd<j which 
enters Eq. (14) and determines dl/dV through Eq. (11). Here we present results using 
three different techniques: a) perturbation theory up to second order in U, b) slave-boson 
mean-field approximation (SBMFA) and c) exact diagonalization plus embedding (EDE). 
The first one has been already used by us to study the mirage effect [38,39,52,53] and by 
one of us [39] and Shimada et al. [80] to study the hne shape of dl/dV in absence of the 
corral. The latter problem was also studied recently using the SBMFA [56] , and to the best 
of our knowledge the results presented in section 8 are the first application of this technique 
to the mirage effect. EDE has been used in Ref. [47]. 



15 



A. Perturbation theory in the Coulomb repulsion 



The starting point is the calculation of the non- interacting problem {U — 0) but with 
Ed replaced by effective one-particle d level E^^-^ . Using equations of motion similar to Eqs. 
(12), and assuming V^G'^lRi, Ri, uj) — —iSb independent of a; (as in the simple case of section 
3), the resulting non- interacting impurity Green's function becomes: 

GdaH = r, • (15) 

uj-E'// + t6,-{VsXrG^,{Ri,Ri,uj) 

The first choice for E^l^ would be the Hartree-Fock value E^^ — Ed+U {d\da) [61]. However, 
out of the symmetric case Ed + U/2 = ep, better results are obtained if E^-^ and {d)„da) are 
calculated self-consistently using interpolative schemes that reproduce correctly the physics 
not only for small U but also for infinite U [23,70,79,81]. For example, the persistent current 
in small rings with quantum dots practically coincides with exact results for U ~ 6A, where 
A is the resonant level width [79] . At the symmetric case the theory is quantitatively correct 
up to [/ ~ 8 A [82]. To avoid self-consistency we take parameters near the symmetric case, 
for which E^JJ is near the Fermi energy ep- This is consistent with first-principle calculations 
[83]. 



< 




> 





FIG. 2. Fcynman diagram for the contribution to the self energy of the d electrons in second 
order in the Coulomb repulsion U. 
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The interacting impurity Green's function can be written in the form 

G^i^) = [G'MV + ^t' - Eaii - S.M, (16) 

and the approximation consists in calculating I^a{^) in second-order perturbation theory in 
U [60,61] (the first-order terms are already included in E^JJ^). The corresponding Feynman 
diagram is shown in Fig. 2. Using the analytical extension of the time ordered G^^^{uj) to 
Matsubara frequencies, the expression for the self energy reads 

E(,(ia;„, T) = U'^T ^ G^^{iuJn - ii^m)x{i^m); 

m 

X{ii^m) = -TY,G%{iuJi)G%{iuJi + IVm), (17) 
I 

where the a;, (z/^) are fermionic (bosonic) frequencies. The evaluation of the Matsubara 
sums is greatly facilitated by the fact that the unperturbed Green's function for surface 
states G°(i?i, i?2,<^), which for a soft corral involve a continuous distribution of energy, can 
be well approximated by a sum over a finite number of simple fractions with simple poles in 
the complex plane [see Eq. (32) of section 6] [52] . 
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FIG. 3. (a) Self energy and (b) impurity spectral density as a function of energy for a system 
with an impurity placed at position (—0.4a, 0) inside an elliptical quantum corral with eccentricity 
1/2 and size such that the state 35 is at the Fermi level. Parameters are in the text (section 5 A). 



In Fig. 3 (a) we show the resulting at zero temperature for an elliptical corral 
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with eccentricity e = 1/2, like that used in the experiment of Manoharan et al. [12], but 
with semimajor axis reduced to a — 6.46 nm so that the state 35 of the hard-wall corral 
falls at the Fermi level. The impurity was placed at a maximum of the wave function of this 
state {x — ±0.4a, y — 0, see Fig. 10). For simphcity we took Vb — 0, and (^^(-Ri, -R2,<^) was 
constructed from hard-wall eigenstates broadened by an imaginary part 5 — 20 meV [Eq. 
(32) ]. For the hybridization with surface states we took an energy dependent decreasing 
function Vs{e) — 0.67eVmax (l -|- ^f^,0) [38], which leads to a more symmetric impurity 
density of states pda{^)- The zero of energy is set at ep — 0, and E'^JJ = 22 meV. We took 
U = 1 eV. While U = 2.84 eV has been estimated [67], this approximation ceases to be 
reliable for larger values of U [84]. The imaginary part of ^^.(a;) vanishes at ep and has a 
quadratic dependence on energy near ep , respecting Fermi liquid properties [85]. 

The particular structure of Eo-(a;) near ep leads to the development of the Kondo peak in 
the impurity spectral density Pda{^^) — — :^IniG^ci(T(<^)- This function is shown in Fig. 3 (b) 
for a range of energies extending between the bottom of the surface band and the smooth 
cutoff in the hybridization. The overall structure is similar to that shown in Fig. 1 (a), 
with two charge fluctuation peaks and the Kondo peak. However, the uneven structure of 
the conduction band, which in this case is a sum of broadened peaks rather than a fiat 
band, introduces some wiggles. This is particularly clear for the charge fluctuation peak 
near E^ + U 0.5 eV. The effects of the confinement will be discussed in the next section. 

Unless otherwise indicated, the results presented in this paper were obtained by this 
technique. 

B. The slave-boson mean field approximation (SBMFA) 

This approximation for the [/ — > oo limit of the Anderson model is in some sense a 
complement of the previous one, which is valid for small or moderate U [84]. The slave- 
boson representation consists in writing dl — f^b as a product of a fermion operator and 
a bosonic one b [66]. For [/ — > oo, double occupancy is forbidden and this is expressed by 
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the constraint b^b + /J-Zo- — 1 introduced by a Lagrange multiplier A in the Hamiltonian 
H [66,86]. We present the formalism for the SU{N) generalization of our model, in which 
the index a can run over a set of N degenerate states (instead of only 2). In mean field, the 
bosonic operators are replaced by a number 6 — > (6) = (b^) — r, and r and A are obtained 
minimizing the free energy of the resulting model for free fermions. In this approximation, 
the charge fluctuation peaks (at and Ed + U) are absent in the spectral density. However, 
in the Kondo regime, for zero or small temperature and energies near the Fermi energy, the 
approximation seems to be reliable [6]. We restrict our calculations to T = 0. 

In the SBMFA, the impurity Green's function near the Fermi energy ei;' = is just 

G^{uj) = r^Gf^{uj), (18) 

and the Green's function Gf^iuj) = {{fa',fa))u) is obtained solving the following effective 
Hamiltonian, which results from Eq. (1) with the above explained replacements 



JO- ja 



+r 



J2^Vmsj, + H.c.) + E(n'/>.> + H.c.) 



Minimization of the energy {H^ff) with respect to A leads to 



(19) 



r2 = l-E^- = l-^^c^> (20) 

(7 

where in the second equality we assume SU {N) invariance and 

Ua = {flU) = -lj duf{u)lmGf,{uj), (21) 

where /(a;) is the Fermi function. 

Using the Hellmann-Feynman theorem [87], the other equation to be solved self- 
consistently reads 

^^^ = A + 5 + S = 0, (22) 

2r or 
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T.(^sns,. + H.c), 5 = 1 5^(V;V]&,. + H.C.). (23) 

The expectation values entering this equation can be evaluated as integrals over f{u) times 
the imaginary part of Green's functions of the same form of the first member of the second 
and third of the Eqs. (12) [88]. Prom the differences between H (Eq. (1)) and H^ff (Eq. 
(19)) one sees that these equations can be used with da replaced by and a factor r 
multiplying Vj and . Then 



S 



-— Im 

TT 



JO- 



u + tr] — e^. 



--V'XHm 



TT 



J dujf{u;)G%Ri,Ri,u;)Gfa{cj) 



(24) 



where for the last equality we used Eqs. (2) and (13). In a similar way one has 



B = -—V^^lm 

TT 



(25) 



while the / electron Green's function is 



(26) 



uj-Ed-k- {rVsX)^G%Ri, Ri, lu) - {rVb)^Gl{Ri, R^, u;) ' 

Prom the self-consistent solution of Eqs. (20) to (26) we obtain Gfa{uj) and r. The 
differential conductance is then obtained using Eqs. (11), (14) and (18). For the calculations 
shown here, we take N — 2 because otherwise the line shape becomes too asymmetric in 
comparison with experiments for the expected 3d^ configuration of the Co impurity [46] . 

In the absence of the corral, for an impurity on a clean surface, we assume constant 
symmetric density of states as in Eq. (4) : 



A Gg{Ri, Ri,uj) — ps 



In 



Ws + UJ 



Gb{Ri, Ri, ^) — Pb 



In 



Wb + U! 

00 -Wb 



(27) 



The self-consistent equations are rather easy to solve for this case and allows us modify the 
parameters to fit the observed line shape. In Pig. 4 we compare the analytical expression 
used by Knorr et al. [17] to fit the low energy part of dl/dV and our results within the 
SBMPA. The same set of parameters is used in section 8 to study the modifications of the 
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line shape in a circular corral. For the bulk density of states we take pb — 0.145/eV per site 
and spin from its value at the Fermi energy reported by first-principles calculations [89] . 
is determined from the filling of one electron per site 2WbPb = 1. From the effective mass 
m* = 0.38 me, where nie is the electron mass [28,90] and a parabohc dispersion, one gets 
Ps — 0.045/eV per site and spin. From the bottom of the surface band we take Wg — 0.4 eV, 
and we assume for simplicity the same value for the high energy cutoff. As mentioned before, 
the results are only weakly sensitive to the cutoff. Ed — —0.8 eV is taken from Ref. [67]. 
The ratio Vb/Vs is determined by imposing a fixed ratio of the resonant level due to bulk 
{Sb = TTpbV^) or surface {5s = TrpsVg) states: a) db = Sg, b) 6b = 36s. The magnitude of the 
hybridization controls the width of the line shape and is a fitting parameter. The value of p 
in Eq. (7) is fixed in such a way that for a clean surface, near 1/2 of the intensity of dl/dV 
is due to bulk states [17,58,59]. Therefore we took p = 1/3 to compensate the approximate 
ratio Pb/ Ps — 3. Instead q is used as a fitting parameter which controls the asymmetry in 
dl/dV. In addition, for the fit of Fig. 4, we shifted the minimum of dl/dV and used a factor 
that represents the quantity ^A^pt in Eq. (10). From the fitting procedure we obtain a) 
q = 0.04, Vs = 0.895 eV and H = 0.499 eV, b) ? = 0.035, Vs = 0.604 eV and H = 0.583 eV. 



1.30 




CO [meV] 

FIG. 4. Comparison between the analytical expression used in Ref. 17 to fit the low energy 
part of dl/dV and our results within the SBMFA. The parameters are explained in section 5 B. 
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C. Exact diagonalization plus embedding (EDE) 



This method developed for impurity problems [62,63], consists in solving numerically 
by the Lanczos method part of the system Hq which contains a finite number of relevant 
many-body states, and treating a one-body term H' which connects it to the rest of the non 
interacting system H^, by an approximate method. For example, H — Hq+ H'+ can 
describe a quantum wire with an embedded quantum dot modeled by the impurity Anderson 
model in a chain [22,24,64]. In this case Hq contains the dot and a few adjacent sites, and 
H' is the hopping of the extreme sites included in Hq to their nearest neighbors in Hj.. For 
an impurity in a quantum corral, Hq should contain the impurity and a few conduction 
eigenstates of the hard wall corral, which acquire a finite width S due to hopping to the rest 
of the system [47]. As we show in the next section, this width is essential to describe the 
physics. 

The method starts by solving the one-particle Green's functions for H' = 0. Those for 
Hr are known, and those of Hq are calculated using the recursion technique combined with 
the Lanczos method. Off diagonal matrix elements are calculated from diagonal elements 
of hybrid states, involving sum and difference of basis states. This information is gathered 
in a matrix g. For a non- interacting system {U — 0), the Green's function of H, which 
we denote by G, can be calculated from the Dyson equation G = g + gH'G. This is taken 
as an approximation for the interacting system. Obviously the approximation is exact for 
H' — and any value of the interaction, and also in the non-interacting case. 

This approximation should be used with caution and incorrect results can be obtained if 
it is applied outside its range of validity. For the Anderson model, a reasonable criterion is 
that the size of the exactly solved part should be smaller or of the order of the characteristic 
length ^ ~ hvp/Tx mentioned in section 4 [22,24,64]. In practice, even when ^ is ten times 
larger than the size of the system, the resulting value of the impurity spectral density at the 
Fermi energy Pdai^p) practically coincides with the exact value, known from Friedel's sum 
rule [22,64] For much larger ^, the approximation is not valid. For example, the width of 
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the Kondo resonance near the symmetric case Ed-\-U /2 ep behaves as U~'^ [64], what is 
incorrect for large U (implying small Tk and large ^) [65]. In the mirage experiment, using 
the velocity of bulk states vp = 1-57 x 10® cm/s, then ^ ~ 200 nm and the size of the ellipse 
is 2a ~ 14 nm. 

This technique is easier to implement than others for the case of more than one impurity 
in the quantum corral and will be used in section 10. 

VI. THE ROLE OF CONFINEMENT 
A. One-body effects 

In the experiments of the mirage effect in an elliptical corral with eccentricity e = 1/2 and 
semimajor axis a = 7.13 nm, the space dependence of AdI/dV reminds the wave function of 
the state number 42 for a two-dimensional free electron gas in a hard-wall corral [12]. This 
already points out the importance of the confinement in the problem. Although the hard 
wall is not a realistic assumption, some basic features of the mirage effect can be understood 
with it [32,33]. The eigenstates which determine the surface Green's function Eq. (13) have 
in general a continuous distribution in energy, but the spectrum is discrete for a hard-wall 
corral. From the form of the Schrodinger equation we know that for corrals of the same 
size, the separation between any two energy levels is inversely proportional to the area of 
the corral [52]. Therefore in principle changing the size of the corral allows to single out 
one energy level at will, place it near the Fermi energy ep, and observe it by STS, since as 
explained in section 4, it essentially captures the conduction states near ep- While this is a 
good starting point for the understanding of the phenomenon, due to the soft character of 
the walls, the corral eigenstates become resonances and there is a delicate interplay between 
the width of these resonances and the separation between energy levels. 

While in presence of soft walls the surface eigenstates form a continuum, it turns out very 
useful not only for the understanding of the physics but for the practical implementation 
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of the many-body techniques, that under general assumptions, the surface Green's function 
can be written as a discrete sum of contributions from resonances [51]. We have shown this 
exphcitly for the case of a circular confining potential of the form 



2m*r5 ro 



V{r,9) = V^^,^^-^d{--l), (28) 



where r, 9 are the polar coordinates on the plane and Konf is a dimensionless constant 
controlling the strength of the confinement potential [52,53]. The result is 

G',{r, e, r', e', a;) = y " Jm[i^nr)Jm[Knn e 

n,m "• " 

The complex poles — iS^^ — (/j.A;^)^/(2m*), where the complex wave vectors /c^ are the 
zeros of a function Fm{k) explained below which lead to positive 5™. The coefficients C™ are 

i k"^ 

dk lfc=A;^ 

Fm{k) is a function of the Bessel functions of the first {Jm) and second (1^) kind which is 
related to the normalization of the wave functions in the continuum. Its expression is 

Fm{k) = Alik) + BUkl with A^{k) = 1 + yJZ^^'tl,,^,, , 

Ym{kro) Jm(kro) 

In practice, in Eq.(29) one includes only the terms for which < Eci where the cutoff 
energy Ec is of the order of 1 eV. 

The results of Ref. [52] suggest that for corrals of other shapes and not too weak con- 
finements, one can approximate the surface Green's function as 

GliRM^^^-^^, (32) 

^—^ LO — €j + idj 

where <Pj{R) are the discrete eigenstates of the hard wall corral and ej are their energies, 
calculated with a shghtly renormalized effective mass m* (increased by about 10% [52]), and 
Sj are the width of the resonances, which to a very good approximation are linear in energy 
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Sj — Spiej — eb)/{eF — Cfc), where 5f is the width of the resonance at the Fermi level and 
is the bottom of the surface band [52] . As we show below, this width plays an essential role 
in the many-body results. In some cases, for simplicity and since it does not affect much 
the results, we will take constant 5j — Sp- 

B. Many-body effects 

Usually, as in the simple case of section 3, the Anderson impurity is hybridized with a 
continuous band of conduction states, flat on the scale of T^. However, if one takes Vb — 
and a hard-wall assumption for the surface states, the Anderson impurity in our model is 
mixed with a discrete set of conduction states with a significant separation between adjacent 
levels. Does a Kondo resonance form in this case? This question has been addressed in the 
context of mesoscopic systems [49]. We illustrate it with a simple problem of a ring of L 
sites described by a tight binding model with hopping t, in which one particular site has 
on-site energy — —t, a Coulomb repulsion U — 2t and hopping OAt with their nearest 
neighbors [79]. We take half fiUing, what implies ep — and assume that the ring is 
threaded by half a flux quantum in order to have (as in the case of the quantum mirage) 
an important hybridization of a conduction state at the Fermi energy ep with the impurity. 
This is a symmetric Anderson model with a discrete spectrum of conduction states. The 
impurity spectral density Pdai'^j) calculated with perturbation theory in U is shown in Fig. 
5. For L — 800, the average separation between the levels which hybridize with the impurity 
d = 8t/L is an order of magnitude smaller than the half width of the resonance Tk ~ O.lt 
and we can see a structure similar to Fig. 1 (a). In particular, the Kondo resonance at 
can be visualized. For L = 80 one has d ~ Tk and the spectral function has some similarities 
with that of the continuous conduction band, but with an important internal structure. For 
L = 8, for which d ~ lOT^:, the central Kondo peak is absent [91]. 
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FIG. 5. Impurity spectral density as a function of energy for a ring of L sites described by 
the impurity Anderson model and several values of L. An artificial broadening 77 = O.Oli was 
introduced for each peak. The result for L = 800 is also shown with a dashed line in (a) and (b) 
to facilitate comparison. Parameters in the text (section 6 B). 

As we see, to obtain a well defined Kondo resonance with discrete conduction states, 
it is necessary that d < Tk- Instead, in the mirage experiments one has d ^ Tk- While 
Tk ^ 5 meV [12], the average distance between the energy levels that have an important 
hybridization with the impurity (those shown in Fig. 10) is of the order of 100 meV. This 
shows the need to take into account the finite width dj to the conduction states. The 
evolution with 5j (taken for simplicity independent of j) of Pda{^) ^-nd the change in the 
surface part of the conduction density of states after addition of the impurity Apsa-{uj) for the 
elliptical corral with e = 1/2 studied experimentally [12] is shown in Fig. 6. The impurity is 
placed at one focus of the eUipse. Psa{^) is given by Eqs. (7) , (11) and (14) with p — q — 0. 
As anticipated above, for very small S, the impurity spectral density does not show a well 
defined resonance at ep. As a consequence, there is a marked disagreement of Apsa{u!) with 
the observed AdI/dV (which is very similar to the bottom left curve). For S — 1 meV, 
Psct(u;) has two peaks (instead of antiresonances as in section 3) at the same positions of 
Pda{^), while in the absence of the impurity Psai'^) has a peak which corresponds to the 
state 42 which hes at ep- Therefore the depression of Apsa{uj) at ep is a consequence of the 
subtraction and does not indicate a Fano antiresonance. These results are consistent with 
numerical results which correspond to 5 = but include (as usual in these calculations) 
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an artificial broadening of the resulting peaks [48]. For a large broadening the results for 
^Psaioj) look like those of Fig. 6 for 5 = 10 meV but with a large positive average, what is 
inconsistent with experiment. 
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FIG. 6. Impurity spectral density Pda{^) (top) and change in the surface spectral density due 
to addition of the impurity at the impurity site Apsa{uj) (bottom) for the configuration of the 
mirage experiment and different values of the width of the conduction levels S : 1 meV (left), 10 
meV (middle) and 50 meV (right). The dashed line shows Apsai^o) at the empty focus. Parameters 
as in Fig. 4 (section 5 A). 

As 5 increases, the two peaks in pda{'^) merge into one (for 5 ~ 18 meV) and the shape of 
both, the Kondo resonance and the Fano antiresonance, become similar to the results of the 
more conventional case, described qualitatively by the simple model of section 3. The Fano 
dip in /S.pstj{uj) for 5 = 50 meV agrees well with experiment [12]. The rather symmetrical 
shape is due to the fact that Vs decreasing with energy was assumed. For constant Vg, 
Apso.(a;) is smaller for positive uj (see full line of Fig. 16). The evolution of poail^) with 
5, described first in Ref. [38], has been confirmed by exact diagonalization plus embedding 
[47] , and by Wilson renormalization group calculations [44] . 
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In Fig. 6 we also show Apsct(<^) at the empty focus. The comparison with the corre- 
sponding value at the focus where the impurity is located establishes the intensity of the 
mirage effect. For small 6 the "transmission" of the Kondo effect to the empty focus is nearly 
perfect, because the space dependence follows closely the density of the state 42 which has 
maxima at the foci (see Fig. 8). As 5 increases, the intensity of the mirage decreases as a 
consequence of interference effects described in the next section. 
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FIG. 7. (a) Impurity spectral density and (b) Apsa{uj) at the impurity site as a function of 
energy for Sj = Spi^j — ei,)/{eF — efe) with Sp = ^ meV. Other parameters are Vg = 0.45 eV, Sj, = 35 
meV, E^^^ = 7 meV, and U = 1 eV. 

The introduction of moderate hybridization of the impurity with bulk states does not 
affect the need to include a non vanishing S to be able to obtain a reasonable agreement 
with experiment. In Fig. 7 we show again both densities for 5 = 1 meV and parameters 
such that the strength of the hybridization of the impurity with bulk and surface states is 
approximately the same [92]. The peaks in padl^) and psail^) are broadened with respect to 
the previous case, but again the dip in /S.psa{oj) is not a Fano antiresonance, but correspond 
to minus the peak in Psaij^) at ej? in absence of the impurity. 
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VII. THE SPACE DEPENDENCE OF DI/DV 



While scattering theories based on a phenomenological phase shift for the scattering at 
the atoms of the boundary and the impurity describe quantitatively the space dependence 
[31,36,37], approaches based on wave functions of a corral (with continuous boundaries) 
usually bring more insight into the underlying physics [32,33,38]. For example the predic- 
tion of mirages out of the foci of elliptical corrals are somewhat hidden in the scattering 
approaches. Instead, mirages observed in a circular corral [13] were inspired by the extrema 
of the wave functions of the degenerate 37*^^ and 38*^^ conduction eigenstates of a hard-wall 
circular corral, and were calculated with our many-body approach for a circular corral with 
soft walls [52,53]. 




FIG. 8. Contour plot of the wave functions of an elliptical corral with semimajor axis a = 7.13 
nm and eccentricity e = 1/2 which lie close to the Fermi energy. 
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FIG. 9. Contour plot of dl/dV in the elliptical corral with a = 7.13 nm and e = 1/2 for S = 50 

meV (left), AdI/dV for S = 50 meV (middle), and AdI/dV for S = 20 meV (right). The apphed 
voltage is 10 mV. Other parameters as in Fig. 4 (section 5 A). 

Having in mind the most studied case of the mirage effect: a Co impurity placed at 
one focus of an eUiptical corral with e — 1/2 built on a Cu(lll) surface [12], we have 
calculated the differential conductance dl /dV{Rt, Ri,V) as a function of the tip position 
Rt, for the impurity position fixed at the left focus [Ri = (—0.5a, 0)] and voltage = 10 
mV. We used Eqs. (11), (14) and (32). We have taken p — q — since they are important 
only near the impurity [74]. Therefore the results depend on the impurity Green's function 
Gd{uj) calculated as is section 4 A, and mainly on the unperturbed surface Green's function 
G°(i?i, i?2, co). To calculate the latter we used Eq. (32) with the corral wave functions <Pj{R) 
calculated as in Ref. [93]. The wave functions of the states which he nearer to the Fermi 
energy are shown in Fig. 8. The wave functions can be classified by symmetry into the four 
irreducible representations of the point group 6*2^, according to the parity under reflection 
through the major (minor) axis ay (ax). In particular each of the shown wave functions 
belongs to a different representation. The state 42, which lies at the Fermi energy is even 
under both reflections, 40 is odd under both of them, 41 is even under ay and odd under 
ax, and 43 is odd under ay and even under ax- 

The results presented in Fig. 9 were obtained for Vb — 0, but quite similar results come 
out if Vs is decreased by a factor 1 / and Vb is increased so that the contribution to the 
resonant level width of bulk and surface states has the same magnitude, and the width of the 
impurity spectral density is kept. This is not surprising since the above mentioned change 
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of parameters practically does not affect Gd{oj)-, and then, from Eq. (14), the only change 
in AdI/dV{Rt) for p = g = comes from a factor (see Fig. 3 of Ref. [52]). Instead, 
the dependence on the impurity position Ri should be affected by the relative strength of 
Vb and Vs (see next section). 

The differential conductance dI/dV{Rt) for a constant width 5 — 50 meV of the conduc- 
tion surface states is represented in Fig. 9 left. It is very similar to the observed topograph 
[12]. However, the latter corresponds to the total current / and not to dl/dV. The similar- 
ity is due to the fact that 5 is larger than the energy corresponding to the applied voltage 
eV = 10 meV, and dl/dV does not change too much in this energy scale. Comparing with 
Fig. 8, one sees that as a first approximation, the observed pattern can be described as a 
sum of the densities of the state 42 which hes at the Fermi energy ep, and the state 43 which 
is ~ eV above ep- The wave function of the state 42 shows some vertical "stripes" which end 
in "arcs" at the extreme left and right. These essential features rotated 90 degrees describe 
roughly the wave function of the state 43. Therefore the structure with "arcs" at the border 
and "stripes" in the middle is to be expected in the sum of probability densities. 

Translated into equations, this is consistent with the behavior expected from the first 
term of the first Eq. (14) and Eq. (32). However, this is not the whole story because the 
above mentioned term depends on the sum of squares of wave functions and is therefore 
invariant under all symmetry operations of the ellipse, while the observed topograph (and 
our calculated dl/dV) has not a defined parity under Ux- The addition of the impurity 
breaks the symmetry under and the effect of the impurity is contained completely in 
AG^j [see Eq. (14)]. The imaginary part of AG/j is directly proportional to AdI/dV, which 
is dl/dV minus the corresponding quantity for the empty corral. From Eqs. (14) and Eq. 
(32), it is clear that physically the effect of this subtraction is to eliminate the contribution of 
all states j which have a negligible hybridization with the impurity V^(^^(i?j). In particular, 
states odd under Oy like 40 and 43 have (p'j{Ri) — and do not participate in AdI/dV. In 
practice, states like 41 which arc even under ay but have a small amplitude at the foci do 
not affect the result either. In Fig. 9 middle we show the "cleaned" result AdI/dV for the 
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same parameters of the complete result dl/dV shown at the left. Now the main features of 
the wave function of the state 42 can be recognized directly, particularly if the width of the 
conduction levels is reduced to 5 = 20 meV (Fig. 9 right). Comparison with experiment 
[12] indicates that the right value of 6 is in between those shown: 20 meV < 6 < 50 meV. 
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FIG. 10. Contour plot of the wave functions of an elliptical corral with semimajor axis a = 7.13 
nm and e = 1/2 with appreciable amplitude and the foci and close to the Fermi energy. 

An analysis of the magnitude of the wave functions at the foci (which determine the 
hybridization strength of the impurity with the different states) shows that the space de- 
pendence of AdI/dV is dominated by four states. The wave functions of these states are 
shown in Fig. 10. While all these states are even under ay (otherwise they would not 
hybridize with the impurity), only 42 is even under ax- The rest are odd under a^- This 
produces a negative interference between the contribution of the state 42 and the other three 
at the empty focus which tends to destroy the mirage effect. In simple terms, one could 
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say that the information of the Kondo effect transmitted by the focus of the impurity by 
the four wave functions reaches the other focus with positive sign for the states 42 and with 
negative sign by the states 32, 35 and 51 so that the amphtude is reduced. Formally, this 
can be seen in Eqs. (14) and (32). As S decreases, the relative contribution of the state 42 
which lies at the Fermi energy increases and the size of the mirage effect also increases. This 
suggest to reduce S, if one can control this parameter experimentally, or to try to optimize 
the geometry in order to reduce the negative interference effects [38]. However, as shown in 
the previous section, if 5 is reduced too much, the Kondo resonance and Fano antiresonance 
near the Fermi level are destroyed. 



VIII. VOLTAGE DEPENDENCE OF ADI/DV 

The experimental study of the line shape of AdI/dV at the impurity site in different 
positions of one corral or in different corrals and its comparison with theory should be useful 
to elucidate the relative strength of the hybridization of the magnetic impurity with bulk and 
surface states. A greater sensitivity to geometry points towards a greater relevance of surface 
states. Recently, it has been argued that due to the exponential dependence of the Kondo 
temperature on the density of states [65], the observed line shape with approximately the 
same width in different situations, indicates that the hybridization with bulk states should 
be much more important [45] . However, the calculations of Ref . [45] are rather generic and 
the specific features of the corral states were not taken into account. 
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FIG. 11. (a) Analytical function used to fit dl/dV for a Co atom on a clean Cu(lll) surface 
(Ref. 17). (b) Density of the state proportional to AdI/dV [see Eq. (11)] for the system of (a), 
(c) same as (b) for a Co atom at the focus of the elliptical corral of Fig. 10, at the Co position 
(full line) and at the empty focus (dashed line). For comparison, the dotted line is the result for 
the clean surface with q = 0. Parameters are Vg = 0.64 eV, S = AO meV, U = 1 eV and q = 0.03. 

For the case of the Cu(lll) surface, to the best of our knowledge the dependence of 
AdI/dV on bias voltage has been reported only in two cases: the clean surface [12,17] 
(see Fig. 4) and the elliptical corral described in the previous section, with a Co atom at 
one of the foci [12]. In the latter case the hue shape is more symmetric, but the width 
is approximately the same in both cases. Both line shapes can be qualitatively described 
including only hybridization with surface states. In Fig. 11 we show our results obtained 
within perturbation theory, using Eqs. (14), with p — 0, and q — 0.03, to control the 
asymmetry of the line. Since we used here constant Vg and surface density of states Ps (as 
in section 3), and the nearly symmetric case Ed + U/2 ~ ep, the line shape for the clean 
surface is symmetric for g = (dotted line in Fig. 11 (c)), and a value of g > reproduces 
the observed asymmetry (Fig. 11 (a)). Instead, a constant Vg in the corral case leads to an 
asymmetry opposite to that observed for the clean surface (like the full line of Fig. 16), while 
the hne shape observed in the corral is symmetric. Then, in this case the effect of g > is 
to correct the asymmetry. It is encouraging that the same set of reasonable parameters can 
explain qualitatively both line shapes. The experimental AdI/dV has kinks around ±0.01 
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V which are probably due to pecuharities of the non-interacting band structure and are out 
of the scope of our theory [12,17]. As shown at the bottom of Fig. 6, the width of AdI/dV 
has some variation with the width of the conduction states 5. Here we have chosen 5 = 40 
meV, which as discussed in the previous section, leads to a space dependence of dl/dV in 
agreement with experiment. 
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FIG. 12. Apft(,(w) - /^dI/dV{u/e) [see Eq. (11)] for an ellipse with a = 6.46 nm and e = 1/2, 
and two impurity positions: (a) Ri at one focus, (b) Ri = (±0.4a, 0). Full (dashed) line corresponds 
to the tip at Ri {—Ri). Parameters are 5j = Spi^j — eb)/{^F — £&) with Sp = 40 meV, Vg = 0.48 
eV, Sb = 32 meV, E^^^ = 7 meV, U = 1 eV, p = and q = 0.02. 

The rather similar line widths in the above mentioned cases seems accidental and other 
situations are more suitable to analyze the relative role of surface and bulk states in the 
formation of the Kondo resonance. In Fig. 12 we show the line shape expected in a smaller 
elliptical corral, with semimajor axis reduced to a = 6.46 nm keeping the same eccentricity 
e = 1/2, so that the state 35 (see Fig. 10) falls at the Fermi energy. This state has extrema 
at positions Re = (±0.4a,0), and the average separation of the levels is larger than in the 
previous case. The surface spectral density near the Fermi level is larger at Rg than at the 
foci (±0.5a, 0). Even including the same hybridization strength of the impurity with surface 
and bulk states, the depth and width of AdI/dV is substantially larger if the impurity is 
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placed at Re rather than at the foci. Note also that the intensity with the tip at the opposite 
point Rt = —Ri is considerable larger in this case. This is due to the fact that the negative 
interference between states 42 and 35 explained in the previous case is substantially reduced. 
A stronger mirage in this geometry has been predicted before [38] . 

In the rest of this section, we show results for the line shape for the tip placed on the 
impurity and several positions of the impurity inside a circular corral of radius tq = 6.35 
nm, in such a way that the degenerate states 37 and 38 lie at the Fermi level. Experiments 
in this corral have been done to illustrate the simultaneous presence of two mirages [13]. 
We use the SBMFA, because it gives the correct exponential dependence of Tk with the 
couphng constant [65]. We also use the exact Eqs. (29) to (31) for the surface Green's 
function instead of the approximate one Eq. (32). The SBMFA is described in section 5 B 
and the parameters are those taken there to fit the line shape for the clean surface. 
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FIG. 13. Apfi^{io) ^ AdI / dV {lo / e) within SBMFA for a circular corral of radius ro = 6.35 
nm at the impurity position {Rt = Ri), for Sp = 40 meV, Sb = Sg and different values of \Ri\/ro 
indicated inside each figure. Other parameters as in Fig. 4 (section 5 B). 

In absence of the impurity, the density of surface conduction electrons at the Fermi 
energy has a pronounced relative maximum near r = \Ri\ = 0.1 5ro [52]. As shown in Fig. 
13 (a), the depth and width of AdI/dV varies considerably as the impurity and STM tip 
are moved together from the center of the corral to this maximum. At the center, only 
the corral surface states with angular momentum projection m = can hybridize with the 
impurity. Since the corresponding resonances are far from the Fermi level (see Fig. 14), the 
Fano antiresonance for r = is more than 80% due to bulk states. In fact, doing the same 
calculation with p = (assuming no hopping between tip and bulk states, see Eqs. (7) and 
(14)), \Apha\ < 0.01/eV. Therefore for r = the bulk states play a major role not only 
in the formation of the Kondo state but also in the variation of the STM current AdI /dV 
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which is mainly due to the current between tip and bulk states. At remote positions this 
Fano antiresonance of bulk states will not be observed [74] . In general, the contribution to 
the dip in AdI/dV due to bulk states (and interference with surface states), captured at the 
impurity by the hybridization of tip and bulk states, will be absent at a mirage point and is 
a natural limitation of the intensity at the mirage point (see next section). For r/rg — 0.15, 
the intensity of AdI/dV decreases to ~ 40% if the tip-bulk hopping is disconnected. 

Compared with the rapid variation for r/ro < 0.2, the width and magnitude of AdI/dV 
oscillates weakly with position for 0.2 < r/ro < 0.9, with larger intensity and width for 
r/ro = 0.4, 0.65, and 0.85 (see Fig. 13). However, there is a dramatic increase for r/ro > 0.9, 
with a maximum near 0.96, as shown in Fig. 13 (e). Although unfortunately at this 
short distances from the boundary our theory ceases to be reliable (because of our simple 
assumption of a continuous boundary potential), it is instructive to relate this result with 
the variation of the density of surface states p°^(r, a;) at the Fermi energy = with 
position. As shown in Fig. 14, there is a moderate increase in p°^(r, a;) near a; = as r/ro 
increases from 0.15 to 0.95. This is mainly due to the contribution of resonances with high 
angular momentum projection which render p^^ rather flat in energy. Instead, for r/ro = 0, 
the resonances with m — are selected and they lead to the displayed oscillatory behavior. 




co[eV] 



FIG. 14. Density of unperturbed surface states p°^(r, a;) = —X'^lm[G^{r,9,r,6,u})/Tr] [see Eqs. 
(29) to (31)] as a function of energy oj for different values oi r/ro. 
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The parameters of Fig. 13 correspond to an equal participation of bulk and surface states 
in the resonant level 5b — Sg. Considering the case 5b — 35^, as expected, the variation of the 
width of the resonance with the position of the impurity is less pronounced, but otherwise 
the same qualitative features as before are obtained. Except for the peculiar behavior near 
the boundary of the corral, the greater sensitivity to the position is for < r/ro < 0.15 
as before. Comparison between both cases is presented in Fig. 15. In Fig. 15 (a) and (b) 
we have used an intensity of the boundary potential Konf = 7, which leads to a broadening 
5f — 4:0 meV of the surface conduction states at the Fermi level (see section 6 A). This 
value leads to a space dependence in agreement with experiment (see section 7). In Fig. 
15 (c) we show how the space variation is affected if the confining potential is increased 
to Vconi ~ 15, leading to 5f — 20 meV. The oscillations in the surface density of states 
and therefore the variation of the width of the resonance with position becomes much more 
pronounced. There is a tendency to a change in the hne shape, similar to that of Fig. 6 for 
5f — 10 meV, which can allow to identify or rule out this regime experimentally. For other 
positions of the impurity (not shown), the tendency is similar to that of Fig. 13, but the 
particular structure for r/ro — 0.95 almost disappeared. 
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FIG. 15. Same as Fig. 13 for (a) Sf = 40 meV, 6b = 36^, (b) 6f = 40 meV, 6b = 6s and (c) 
dp = 20 meV, 6b = 6s. 



IX. LOWER BOUND FOR IMPURITY-SURFACE HYBRIDIZATION 

Within our local picture for the hybridization of the impurity and tip with conduction 
states, a simple estimate of a lower bound for Vg/Vb can be obtained from the mirage 
experiment in the elliptical quantum corral [12]. The ratio of the intensity of AdI/dV at 
the mirage point (for the tip position Rf — —Ri) to that at the impurity /j (for Rt — Ri) 
was reported to be Im/h — 1/8. As in section 3, let us approximate G^{Ri,Ri,uj) ~ 
—mpb- Also, for enough broadening of the surface conduction states 5 > 40 meV one has 
G^{Ri, RijCu) ~ —i'Kps [38]. Neglecting the tip-impurity hopping (g = 0), using Eqs. (11) 
and (14) one has 



h^-MIIdV{Rt)^c{^^+p^^ p^{u)^c{^^+l^ Pdiiv), (33) 

where C = T^VbPs ^ constant and in the last equality we assumed that ppb — Ps, so that in 
the absence of impurity, the tip detects bulk and surface states with the same intensity, as 
reported experimentally [17,58,59]. Now, at the mirage point one can neglect G^(i?i, —Ri^uj) 
[74]. Assuming instead perfect transmission from the surface states (as if only the state 42 
were relevant), one has G°(-Ri, —Ri,uj) = G^g{Ri, Ri,uj). Since the amplitude at the mirage 
point is less than that for perfect transmission, one has from Eqs. (11) and (14) 



Im = -AdI/dV{-Rt)<C\^-^J PaH, (34) 
and then, from Eqs. (33) and (34) 

Solving a quadratic equation, a more precise bound for Im/h = 1/8 gives Vg/Vb > 0.547. 
Using pb — 3ps, this implies 6s > Sb/lO, with Sc = TrpcV^. A smaller tip-bulk hopping p 
leads to a larger lower bound for Ss- 
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X. INTERACTION BETWEEN KONDO IMPURITIES IN A QUANTUM 

CORRAL 

Experiments with two impurities inside an elliptical quantum corral have been done [13], 
but the results were not published yet. These experiments should be particularly useful as 
a test of the relative strength of the hybridization of the impurity with bulk and surface 
states, since one expects that at distances larger than ~0.5 nm, the interaction between 
two Kondo impurities is dominated by surface states. In this section we extend previous 
calculations of the line shape of AdI/dV when there is one impurity at each focus of an 
elliptical corral, using the technique of exact diagonalization plus embedding, described in 
section 5 C [47]. Other calculations of the interaction between magnetic impurities in a 
corral have been made by perturbation theory in the Kondo coupling [94]. However, this 
technique does not work in the case we are interested of antiferromagnetic Kondo coupling 



As explained in section 5 C, the Hamiltonian is written ets H — Hq+ H'+ H^. In 
our case the Hilbert space of Hq contains one or two impurities and the most important 
surface conduction states (those represented in Fig. 10) and two additional ones (24 and 62) 
although they do not affect the results. Hj. describes a set of independent non-interacting 
bulk states which hybridize independently with the impurities and the surface conduction 
states. Hq has a similar form to Eq. (1) but only contains hard-wall surface conduction 
states and can contain more than one impurity: 



Hq^Y, ^A^3- + EdJ2 4adia + UY, 4^di^4idii + J2 ^Vs[<Pj{Ri)dlsj, + H.C.], (36) 



while the effect of bulk states and the broadening of the surface states (necessary to obtain 
a qualitatively reasonable line shape as shown in section 6) is contained in H' which reads 



For the bulk states bi„ we take a constant unperturbed density of = 0.05 states/eV per 
spin (of the order of the density of bulk s and p states at e^. [90]), but a change in pf, can 



Jk>0 [65]. 




(37) 
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be absorbed in a change in t and Vb. The value of t controls the width of the conduction 
states and therefore, the intensity at the mirage point, as explained in sections 6 and 7. 




1 imp 

1 imp mirage 

2 imp 

2 imp V=0 
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FIG. 16. AdI/dV as a function of applied voltage calculated with the EDE for an elliptical 
corral with one or two impurities at the foci. Full line: one impurity and tip on it. Dashed line: 
one impurity and tip at the other foci. Solid squares: two impurities and tip on one of them. Open 
circles: same as before with Vs increased by a factor \/2 and V;, = 0. Unless otherwise indicated 
parameters are Ed = -0.8 eV, C/ = 3 eV (from Ref. 67), Vs = 1. 12/^2 eV, H = 1.2 eV,p = q = 0. 

In Fig. 16, we represent the change in differential conductance for p = q = for the 
case in which there is one impurity at each focus, comparing two situations. In the first one 
Vb — and 1.12 eV is taken to reproduce approximately the experimental width. In 
the second one Vg is reduced by a factor l/\/2 and Vb increased to 1.2 eV in order that the 
same width and practically the same line shape is obtained. In the former case, the effect of 
the interaction between impurities is stronger and the line is wider and with some structure 
due to a partial splitting of the Kondo resonance. In the second case, the result is very 
similar to the sum of the spectra at both foci when only one impurity is present. This is the 
result expected for weak interactions. 
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XI. SUMMARY AND DISCUSSION 



Using a simple impurity Anderson model in which the particular structure of the surface 
states inside a corral is taken into account appropriately, the basic physics of the mirage 
experiments in quantum corrals [12,13] can be understood. The resulting space and voltage 
dependence of the differential conductance dl/dV is in good agreement with experiment. 
The voltage dependence of dl/dV observed for one Co atom on a clean Cu(lll) surface and 
for a Co atom at the focus of an elliptical corral built on that surface can both be explained 
with the same set of parameters (see Fig. 11). 

While the space dependence of dl/dV is mainly determined by non- interacting conduc- 
tion electron Green's functions, the calculation of the dependence of dl/dV on bias voltage 
is a non trivial many-body problem in which the particular structure of the conduction 
electrons for surface states introduces additional complications. Single-particle scattering 
theories were successful in explaining the space dependence of dl/dV [31,36,37], but the 
voltage dependence is actually used to adjust a phenomenological energy dependent phase 
shift. Therefore the differences in the line shape as a function of the apphed voltage in 
different structures (like the above mentioned between an impurity on a clean surface or 
inside the corral) cannot by accounted for. In contrast, the many-body treatment is very 
difficult to implement for open structures, while one-body scattering theory assuming sim- 
ple interactions, allows not only to calculate but also to optimize open structure to obtain 
multiple mirages or other desired effects [95] . On the other hand, our approach leads to very 
good agreement with the space dependence of dl/dV (see Fig. 9) except perhaps for the 
finest details which we did not attempt to fit. In addition, it allows to understand the basic 
observed features, including the mirage effects and its intensity, in terms of the interference 
of wave functions in the corral (see section 7) . 

Experiments in which the change in differential conductance after addition of one im- 
purity in a quantum corral AdI/dV is measured should be able to discern the relative 
importance of surface and bulk states in the formation of the Kondo singlet. Measurements 
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in circular corrals, easier to handle theoretically, would be useful. In section 8 we presented 
some results for this case. In addition, the change of AdI/dV when more than one impurity 
is present in the corral is very sensitive to the surface-impurity hybridization. The obser- 
vation of a quantum mirage establishes a lower bound for this hybridization which we have 
estimated. 

The calculations presented here are for T = 0. Within perturbation theory, the results 
can be extended easily to T 7^ 0. Some results were presented in Ref. [52]. Important 
changes occur in the scale of the Kondo temperature, but the behavior is similar to that 
already known for the simple case explained in section 3. 

An improvement of the many-body theory requires a better knowledge of the hybridiza- 
tion of the impurity with surface and bulk states and their wave vector dependence (which 
we have neglected). The wave vector dependence is expected from a jellium model [55]. In 
particular for the states near the Fermi energy ep and wave functions decaying as exp{—Kz) 
out of the surface, one expects /c^— constant and therefore a smaller wave vector par- 
allel to the surface /cy implies a weaker decay rate k. in the perpendicular z direction, and 
therefore a larger hybridization with impurity and tip. For surface states near ep, is 
small and one would expect a weaker decay for surface than bulk states. This suggests a 
stronger relative hybridization of bulk states with the impurity in comparison with the tip, 
because the former is closer to the surface. This would be consistent with the fact that the 
hybridization of the tip with surface and bulk states is of the same order [17,58,59], and 
the proposal that bulk states dominate the hybridization with the impurity [17,18,45,54]. 
However, more detailed calculations of the matrix elements found oscillations of the matrix 
elements with z, and a stronger relevance of surface states in the formation of the Kondo 
state [56] and in the distance dependence of the observed dl/dV for magnetic impurities on 
clean (111) surfaces [56,57]. 

We have shown that the width of the surface conduction electrons S plays a crucial 
role in the many-body theory (see section 6). We have calculated this width for a circular 
confinement potential and found that it increases linearly with energy [52]. However, we 
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assumed that for a clean surface the surface states are well defined for all wave vectors and 
this is not the case for energies above the Fermi energy [30]. Since localization involves a 
participation of all wave vectors, there is a contribution to the width brought by the states 
of larger wave vector which we have neglected. However, since most of the physics depends 
on the value of S at the Fermi energy and we took it as a parameter, our conclusions are not 
affected. 

Correlation functions of impurities inside a quantum corral have been studied previously 
[47,48,52,94,95]. For perfect confinement a strong enhancement should occur. However, for 
realistic broadening of the surface conduction states, and distances of the order of several 
nm involved in the mirage experiments, we expect that the single ion physics dominate the 
RKKY interactions, and no significant magnetic correlations are present [52] . 
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